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Abstract Three particles floating on a fluid surface define a triangle. The aim of this paper is to 
characterise the shape of the triangle, defined by two of its angles, as the three vertices are subject 
. to a complex or turbulent motion. We consider a simple class of models for this process, involving a 

combination of a random strain of the fluid and Brownian motion of the particles. Following D. G. 
Kendall, we map the space of triangles to a sphere, whose equator corresponds to degenerate triangles 
with colinear vertices, with equilaterals at the poles. We map our model to a diffusion process on 
\ the surface of the sphere and find an exact solution for the shape distribution. Whereas the action of 

the random strain tends to make the shape of the triangles infinitely elongated, in the presence of a 
Brownian diffusion of the vertices, the model has an equilibrium distribution of shapes. We determine 
here exactly this shape distribution in the simple case where the increments of the strain are diffusive. 



Keywords Diffusion, advection 
PACS 05.40.-a,05.45-a 

1 Introduction 



Consider motion of three tracer particles floating on the surface of a fluid which exhibits a complex 
flow JTJ[2] - The three points form a triangle whose shape is specified by two angles, and is therefore 
, specified by a point in a two-dimensional manifold. It is of interest to consider the evolution of the 

shape of this triangle, independent of its scale, rotation and translation of its origin, because this is a 
means to characterise the essential geometry of the flow. Early investigations of random shapes were 
mathematical [3j[4j[5] , and originally motivated by discerning whether nearly co- linear arrangements of 
archelogical features occured by chance or by design [5]. In the context of fluid turbulence, models for 
the evolution of shapes of triangles and tetrahedra have been extensively studied in the past decade 
in the fluid literature: see [TUB ,7,8,9, ID] . Our objective here is to define and analyse the motion of 
triplets of points in a simple model of turbulent flows. We obtain an exact expression for the steady- 
state statistics of the triangle shapes. 
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Our model flow is motivated by considering a scale decomposition of turbulent velocity fields. The 
flow at scales much larger than the size of the triangle merely translates and rotates the three points 
of the triangle, without any shape modification. Structures existing at a scale comparable to triangle 
size are modelled by applying a random matrix of unit determinant, to respect the incomprcssibility 
constraint. We refer to this motion as random strain (although the motion involves both strain and 
vorticity). Lastly, there may also be small-scale components in the flow which result in apparently 
uncorrelated motion of particles. We solve a model incorpoarating these effects, which is a simplification 
of a model introduced in [TJ|S] . Specifically, it is assumed that the particles evolve under the combined 
influence of a random strain (delta-correlated in time) and of independent Brownian motions. Because 
the displacement caused by a uniform strain is proportional to the scale size of the triangle, the strain 
will become dominant in the long-time limit as the separation of the particles grows. We therefore 
analyse a model where the Brownian displacement scales with the radius of gyration of the triangle. 
That is, we analyse a model for the small displacements Sxi of three particles with positions Xi which 
has the following structure: 



where <5A is a traceless random symmetric matrix with diffusive fluctuations, 5ri is a Brownian dis- 
placement, and R is the radius of gyration, defined by 



The stochastic incrcmemts 8 A and Sxi in equation ([TJ are both multiplied by terms which depend 
upon the coordinates X{. This implies a potential source of ambiguity about whether the increment 
5xi at each time step should be evaluated using the Ito or the Stratonovich rule jTTj . We will show 
that this ambiguity only affects the way the size of the triangle grows in time. When considering the 
evolution of the shape of the triangle, we do not need to distinguish whether the stochastic increments 
are of Ito or Stratonovich type. 

In this article, we begin by investigating separately the action of the strain and of the Brownian 
terms. We then consider the two terms together, and find an exact expression for the distribution of 
the shapes of the triangles which are generated by this random process. 

Specifying the shape of a triangle requires two parameters, which could be two of the three angles. 
For our purposes, however, another parametrisation is much more useful. D.G. Kendall argued that it 
is natural to describe the shape of a triangle by a point on the surface of a sphere [3] , which we will refer 
to as the Kendall sphere. His paper also states that if the vertices of a triangle execute independent 
Brownian motions, then the image of this process is uniform diffusion on the surface of the Kendall 
sphere. His original paper [3] does not give a derivation of this result, and derivations which we have 
seen in later works rely upon computer algebra packages [5] or upon quite abstract machinery [3J. In 
section [5] we start by explaining the Kendall sphere, its relation to other parametrisations, and we 
describe a Langevin equation which will prove to be a convenient representation of diffusion on the 
surface of a sphere. In section [3] we define the ensemble of random strain SA introduced in ([!]), and 
evaluate the image of this random strain process on the surface of the Kendall sphere. We show that 
this leads to a diffusion process which is axially symmetric on the sphere. In section|3]we use the same 
framework to give a direct demonstration of Kendall's result that the image of Brownian motion of 
the three points in Cartesian space is unbiased, isotropic and homogeneous diffusion on the Kendall 
sphere. 

In section[S]we combine the results of sections [3] and 0] to determine a Fokker-Planck equation which 
determines the probability density for the model ((T|) on the Kendall sphere. We solve this equation for 
the steady-state shape distribution. 

The analysis in sections[3][S]is uses the strategy of determining stochastic increments of polar coordi- 
nates on the Kendal sphere, and using these to surmise the Fokker-Planck equation. A complementary 
approach is to develop a Fokker-Planck equation for the model ([1]) directly, and make a change of 
variable. This approach is explored in section El where we also show how the Fokker-Planck operator 
can be expressed as a linear combination of Casimir operators. Last, we briefly summarize our results 
in Section [7] 



Sxi = SAxi + RSri 



(1) 




(2) 



i=i j=i 
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2 The Kendall sphere 

Consider a triangle in the plane formed by three points Xj, i = 1,2,3. The shape of the triangle is 
described by two two-component vectors u\ and it 2 and by a 2 x 2 matrix w, defined by 

ui = — t= (xi - x 2 ) , u 2 = — ^(2x 3 - xi - x 2 ) (3) 



and 

' "1,1 "2,1 
"1,2 "2,2 



W 



(4) 



where "i >a , a = 1, 2 are the Cartesian components of Ui. It is natural to represent the matrix w by a 
singular-value decomposition, writing 

w = R( x )diag(A 1 ,A 2 )R(0/2) (5) 

where R(#) is the matrix describing a rotation by an angle 9: 

R(6)=( COse - S ' m °) . (6) 
v ' \ sin cos 1 v ' 

Note that the angle x m © describes an overal rotation of the triangle, and is irrelevant to describing 
its shape. The radius of gyration R of the triangle may be expressed in terms of the singular values Ai 
and A 2 : 

3 

R 2 = \xi - x\ 2 = u\ + u% = tr(w T w) = X 2 + \\ (7) 

i=l 

where x is the centre of mass of the three points. The area of the triangle is simply expressed in terms 
of Ai and A 2 as: 

C = det(w) = AiA 2 (8) 

We note that the size of the triangle can be characterized cither by R or by (. Using C is natural 
when considering area-preserving dynamics [7]. However, because the parametrization involving £ has 
a singularity when the three points defining the triangle become co-linear, it is more natural to use R 
instead of ( when the dynamics allows the vertices to become co-linear. For this reason we prefer R 
over £ as the scale coordinate throughout most of this paper. The angle 4> characterizes the rotation 
in the space of triangles with a given set of values of Ai and A 2 . We note that a cyclic permutation 
of the vertices Xi, (1,2,3) — > (2,3,1) (respectively (1,2,3) — > (3,1,2)) amounts to a transformation 
— » cj) + 27r/3 (respectively 4> — > 4 1 + 47r/3). 

The shape of the triangle may be specified by the variables cj) and 

A =T ^ 

instead of by two angles. 

It is useful to consider the form of the matrix w T w: 

T f U 2 Ui ■ U 2 

w w = I o 

U\ ■ lt 2 U2 




= R(-</ ) /2)diag(A 1 , A 2 )R(- X )R(x)diag(A 1 , A 2 )R(</)/2)) 
= R(-0/2)diag(A?,A 2 ; )R(0/2) 

/1 2 -1 j. \ 

(10) 

so that the geometrical properties of the vectors U\ and u 2 have period 27r in the angle <p. 
It is convenient to transform the ratio A into another variable, z, defined by 

_ 2dct(w) 2A 2C 
2 tr(w T w) 1 + A 2 R 2 ' 1 } 
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Note that the variable z, denoted w in [U[7], satisfies — 1 < z < +1, so that a point the region |z| < 1, 
< <t> < 2ir uniquely defines the shape of a triangle. The points z = ±1 correspond to equilateral 
triangles, with their shape being independent of the angular coordinate <f>. We can therefore regard the 
space of triangles as being equivalent to the surface of a sphere with azimuthal angle <f>, and with polar 
angle 8, related to z by z = cos 8. The equator corresponds to degenerate triangles consisting of three 
colinear points. The poles are equilateral triangles. There are three singular points on the equator, at 
which two points of the co-linear denerate triangle coincide. This device of representing the shape of a 
triangle by a point on the surface of a sphere was first introduced by D. G. Kendall [3J. We will refer 
to it as the Kendall sphere. 

In the case where a triangle is subjected to a sequence of uncorrelated infinitesimal perturbations, 
the representative point diffuses over the surface of the sphere. We will characterise this diffusion 
process by specifying the stochastic increments Sz and Scj) of the coordinates z and <p. In general 
the diffusion process may be inhomogeneous, it may have an anisotropic diffusion coefficient, and it 
may have drift terms. As a reference model, let us consider the case of homogeneous and isotropic 
diffusion on the Kendall sphere, with diffusion coefficient D. We characterise this process in terms of 
the stochastic increments of z and (f>. 

To characterise this homogeneous diffusive evolution let us use a system of standard Cartesian 
coordinates: (x,y, z) = (sin 8 cos 0, sin 8 sin^, cos 8). Without loss of generality we can use invariance 
with respect to the azimuthal angle <p, and consider a point where y = 0. Consider an orthogonal basis 
defined at the point (sin 8,0, cos 8): 

n = (sin 8, 0, cos 8) = (x, 0, z) 
m = (— cos 8, 0, sin 8) = (— z, 0, x) 

k = (0,1,0). (12) 

The diffusive fluctuation of a point n = (x, 0, z) on the surface is 

Sn = SXm + SYk + SZn (13) 

where SX and SY are angular increments, satisfying (SX) = (SY) = and 

(SX 2 ) = (SY 2 ) = 2DSt , (5X5Y) = . (14) 

The value of SZ is determined from SX and SY using the constraint that n remains normalised: this 
implies that SZ = — \{8X 2 + SY 2 ). Consider the increments of z and <f> under this process: we have 
SY = sm85(f> and (n + <5n).e3 = z + Sz = cos 8 — sindSX + cos 8SZ. Solving for S(j>, Sz, we find 

S<P = -SY 

x 

Sz = -xSX - \z{SX 2 + SY 2 ) (15) 



where z = cos 8, x = y/1 — z 2 . The diffusion and drift coefficients for homogeneous diffusion on the 
sphere are defined by writing (Sz) = v z 8t, (Sz 2 ) — 2D z St, (5(f> 2 ) — 2D<f,St, are therefore 

D z = (1 - z 2 )D 
D 

D * = — 

v z = -2zD 

v x = 0. (16) 



This representation of homogeneous, isotropic diffusion on the surface of a sphere will prove to be an 
instructive comparison case for the random strain model which we consider next. 



3 Image of random strain on the Kendall sphere 

Here we consider the dynamics of the three verices fj of the triangle under the action of a matrix M 
which satisfies dct(M) = 1. We shall describe the shape of a triangle by the matrix used to obtain it 
by deformation of an equilateral triangle. Call this matrix M. We consider area-preserving flows, so 
that we require dct(M) = 1. The equilateral triangle (with unit radius of gyration) is described by 

A general triangle is described by 

w = Mw . (18) 
We use a singular- value decomposition representation of M: 

M = R( X ) diag(Ai,A 2 ) R(^/2) (19) 

with A1A2 = 1. 

The strain term, introduced in ([T]), acts multiplicatively on the variables defining the shape M, in 
the form: 

M' = (I + SA) M = R( x + S X ) diag(Ai + SX U A 2 + SX 2 ) R(</>/2 + 5<j>/2) (20) 
where 8 A is an infinitesimal random strain, of the form 

* A =(£-!)• < 2i > 

The statistics of the elements of the random strain are (da) = (Sb) = and 

(da 2 ) = 2DJt , (5b 2 ) = 2D b 5t , (SaSb) = . (22) 

Imposing the reqirement that the statistics of this random strain matrix are rotationally invariant leads 
to the condition D a = Di, = D s . An infinitesimal rotation could have been included in the definition 
of the random strain (by making the off-diagonal elements dissimilar) , but this in clearly irrelevant to 
shape fluctuations. 

Because the noise term 5 A acts multiplicatively on the matrix M defining the shape of the trian- 
gle, the precise definition of <5AM matters when practically integrating the equations. The simplest 
procedure consisting in taking M(i + St) — M(t) « SA ■ M(i) (the Ito procedure), and the physically 
more correct procedure, consisting in defining: M(t + St) — M(t) ~ ^SA ■ (M(t) + M(£ + St)) (the 
Stratonovich procedure) [llj differ by a diagonal term, 2D s 5t I, which dilates in a similar way the two 
eigenvalues Ai and A2, thus leaving the variables that characterize the shape, z and <p, unaffected. For 
this reason, the precise way of defining the model here turns out to be immaterial. 

Consider the evolution of the parameters </>, x, Ai, A2 under this random strain evolution: from (|20l) 
we obtain 

R(S X ) diag(Ai + 5X U A 2 + S\ 2 ) R(5<f>/2) = diag(A 1; A 2 ) + <5A'diag(Ai, A 2 ) (23) 

where 

*A' = R(- X ) SAR( X ) = (% ^ • (24) 

In the following we use the rotational invariance of the statistics of SA to conclude that the statistics 
of its elements satisfy (5a' 2 ) = (5b' 2 ) = 2D s 5t. Expanding the factors of equation ([2"3"]l up to quadratic 
terms in the small increments 5\, S<p and 5\, and collecting terms, we find 



Ai<W 


= 5X 1 -±(46 X 2 


+ Scf) - 




X 2 5a' 


= SX 2 -^(4S X 2 


+ 5<f? ) - 


Y 5 xH 


X 2 5b' 


= A 2 <5x + + 


^XM 


f 5X 2 5x 




= -Ai<5x - 


- SXiSx 


- -5X 2 5. 



(25) 



6 



The first order relations between the increments 8a' , 8b' and the increments of the shape parameters 
are therefore 

8\i = X\8a' 
8X2 = — \28c1' 
4AiA 2 



X 2 \ 2 

2 
2 



8b' 



X 2 + X 2 
\ 2 \ 2 



S X = TT-^5b' (26) 



so that 



(Sx 2 



2D s 8t(X\ + A 2 ,) 2 

{X\- xi) 2 

Z2D s 8tX\Xl 
(A? - Al)2 



(A?-A§) a 

Adding the corrections arising from the quadratic terms then gives 

Ai<5a' = 8X\ — X\ ry— , 2 2 D s 8t 
Aj - A 2 

— A 2 i5a' = 8X2 + A2-T2 TT^sSt 

Aj — A 2 

A 2 £&' = A 2 <5x + ^Ai<50 

A 1( ft>' = -Xx8 X - \x 2 8<l> . (28) 

The first two of these relations are consistent with the constraint <5(AiA 2 ) = 0, and the second pair 
follow from (8a' 8b') = . 

Now transform to new variables A = A1/A2, R = \/X\ + A|- It is clear that the evolution is scale 
invariant, so that the dynamics does not depend upon the radius of gyration R. We find 



6A1 _ A SX2 + ^5A| _ 8X^X2 
A 2 A 2 X\ X\ 



The increment of d> is 



Y A 2 -A 2 . A 2 -l 
The drift and diffusion coefficients in the A, <p variables are: 

D A (A) = 4A 2 D S 
16A 2 D, 

D^(A) = 



2A8a' + 8 f D f 2 t . (29) 
A 2 - A 2 



4A,A 2 w _ AA sy _ 



{A 2 - l) 2 

»0 = . (31) 
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Now work in terms of z = 2 A/ (A 2 + 1): the diffusion and drift coefficients are 



dz d 2 z „ , -MD S A 3 

■■W = ita va + i^ d « a = wtw- (32) 



Then: 



and 



finally 



sin 

A = — , z = sm(20) . (33) 



To express these in terms of z, consider an intermediate variable 0, satisfying 

sin 

COS0 

16D S A 2 {1- A 2 ) 2 



{l + A 2 Y 
AD S sin 2 (2(9) cos 2 (2(9) 

4D s z 2 (l~z 2 ) (34) 



16A 2 D S 4D S sin 2 (26) 

u *\ z )- A 2_ 1 COS 2 (20) 

±D s z 2 
l-z 2 



. -Z2D»A 3 
v z (z) = 



(35) 



(f + A 2 f 

= -8D s z 3 . (36) 

It is interesting to note that the coefficients differ from those of the homogeneous diffusion on a sphere 
(considered in section^ by a simple factor, namely 4z 2 . 



4 Image of Brownian motion on the Kendall sphere 

In the Introduction we mentioned that if the vertices execute independent Brownian motion, then the 
dynamics of the representative point on the surface of the Kendall sphere is homogeneous diffusion. 
The published derivations [U[S| are difficult to compare with the derivation of the diffusion process 
for the random strain model which we analysed in section [31 In this section we present an elementary 
derivation of Kendall's result, using the same approach as for the random strain model. 

Let the vertices of a triangle make small, independent and isotropic diffusive displacements, as if 
they are undergoing Brownian motion. The vectors Ui undergo displacements with components 

5u\ >a = ~^^ Xl - a ~ ^ X '2,a) . SU2, a = -y={26x 3>a - 6x lja - 5x 2: a) ■ (37) 

The statistics of these small displacements satisfy 

(Suij) = , (Sui t j5uk,i) = 25 lk S jl D b 5t (38) 

where is a diffusion coefficient. Note that (8u 2 •) = (Sx 2 t ) for any choice of i, j, k, I, so that is 
the same as the diffusion coefficient for the Brownian motion in coordinate space. 

We can characterise the change in the triangle by considering a change matrix w and hence the 
changes in the variables <fi and z. Applying rotation matrices to w + (5w, we obtain 

R(- x )(w+<Sw)R(-0/2) = diag(A ll A 2 )+R(-x) S £*A R(-0/2) = diag(Ai, A 2 )+<5w' . (39) 
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The elements of 5w' are independent of % and 4>: for example 

Swu = cos x cos(0/2)<5«i,i — sin x cos(0/2)(5wi ;2 + cos x sin(0/2)(5u 2 ,i 



LXsin(0/2)<5 M2 , 2 (40) 



has a variance {{Sw'n) 2 ) = ((Swu) 2 ) which is independent of <fi and x-i an d the same is true of the 
other components of Swij . Because of this invariancc it suffices to perform perturbation theory about 
the x — 0, <f> = case: 



Ai + Su[ j Su' 2 1 

Sib'! J A 2 + Su'r, 



cosSx sini5x 
— sin<5x cos Sx 



\i+S\i 
A 2 + <SA 2 



cos(<ty/2) sin(c50/2)\ 
-sin(4/2) cos(^/2) J 



(41) 

where the Su^ 1 are components of Sui j in the rotated coordinate system. Expanding in powers of 59 
and Sx, and retaining terms only up to quadratic order in the small increments, we find 



<K,i 


= 5\ 1 -±(48 X 2 


+ <5</> 2 ) - 


—5x5(f) 


8u' 22 




+ 54?) - 




5u' l2 


= -Ai<5x - y<50 


(5A 2 (50 
2 


- (5Ai<5x 






5Xi8(h 


- 5X 2 Sx ■ 


Su' 21 


= y 5 ^ + A2 ^ x + 




2 4 



(42) 



These equations can be solved to determine 5Xi, 69 and Sx in terms of the elements of <5w'. We follow 
the same approach as in section [3] First, compute the linear terms from the linear terms of (1421) . We 
then use these expressions to estimate the covariances of the small increments 5<p, 5x and SXi . We find 
that the stochastic increments arc 



SXi = Su'n + 



sx 2 = 



2X 1 D b St 

A l ~ A 2 

2X 2 D b St 



'2,2 



A? 



x 2 

A 2 



2(X 1 Su' 2A + X 2 Su' 12 ) 



Sx 



XiSu[ 2 + X 2 8u' 21 



Ai 



A? 



(43) 



From now on we concentrate upon the evolution of the coordinates which define the shape: these can 
be taken to be <fi and A = Ai/A 2 . Note that 



SA 



Ai + 6X 1 
A 2 + SX 2 



hi 

A 2 



SXi 



A 



6X2 
A 2 



A 



6X1 6X^X2 



Ai 



\2 
A 2 



(44) 



Now normalise the scale of the triangle so that the radius of gyration yAf + Af 
the stochastic increments: 



A 2 



2AD h St 



A? 



Ai 



1 Su[ A Sv' 2 
R A 2 i? 

1 (5^; 



2vl J D, ) <5i 2AD h St 



A 2 
2/H) b <ft 



A? 



i? A 2 i? 
A Sv' 2 



Ai 
2i? 2 



\2 

A 2 

R 2 



A? 



Ai 



R 2 

1 2/1(yT 2 + l) 2 A,<ft 



Ai 



RX 2 R RX 2 R R 2 



A 2 - 1 



R is factored out of 



(45) 



We consider the evolution of the shape of the triangle with the radius of gyration frozen. When the 
radius of gyration is R = 1, the increment of A may therefore be written 



6 A = (1 + A 2 )8x + 



2/1(1 + A 2 ) 2 D b 8t 



(46) 
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where 



5u\ i — A8u'n o 

5x = M 2 ' 2 (47) 



is a stochastic increment which satisfies (Sx 2 ) = 2D b 6t, (6x) = 0. 
Similarly, the increment of <f> may be expressed as 



2VA T +T . 

— T2 — ~°y 



" A 2 (^ 2 -l) 

i 2(A 2 + 1; 

~ i? A 2 - 1 

where 



<fy (48) 



vM 2 + 1 

is a stochastic increment with the same statistics as 5a;, that is (Sy 2 ) = 2D b St. 

The shape fluctuations are described by fluctuations of just two variables, </> and A, which evolve 
diffusively. Their coefficients of diffusion and drift are independent of <j>: 

MA 2 + l) 2 
D A (A) = (A 2 + l) 2 D b 

, AS 2ylM 2 + l) 2 

V4, = . (50) 

An alternative coordinate is z = 2A/(A 2 + 1). Expressing the stochastic increments in the (z,(f>) 
coordinate system, we obtain the following expressions for the diffusion and drift coefficients: 

fdz\ 2 fA 2 -l\ 2 

dz d 2 z _ -16A „ 
V * = dA VA + dA 2 ° A = I 2 TI ' (51) 

Again following the approach in section [31 we express the diffusion and drift coefficients in terms of 
z by introducing a variable defined by A — sin 0/ cos 0, so that z = sin(26>). Now express the 
coefficients in terms of z, and obtain 



_ A lY D _ 4 (sin 2 - cos 2 6>) 2 
iTA 2 ) b ~ (sin 2 + cos 2 0) 2 



D ^ z ) = 4 TT^i D b = -T~^ -T7^T D b = 4cos z (20)D b = 4(1 - z 2 )D, 



fA 2 + l\ 2 16L> 2 4D b 

v z (z) = ^^D b = -8zD b . (52) 



These are the same as for simple diffusion on a sphere, as discussed in section [2J with the diffusion 
coefficient D b replaced by AD b . These coefficients are also very closely related to those obtained in the 
random strain model: the coefficients for that model all differ by a factor of z 2 . 
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5 Mixed model 

Consider the distribution of triangle shapes for a mixed model, in which the displacement of the three 
particles in a short time interval St is a combination of a random strain and an uncorrelated random 
displacement, which is proportional to the radius of gyration R. This leads to the following map, using 
the notation of (U}: 

x'i = Xi + Sxi = (I + SA)xi + R8r l (53) 

describing the evolution of the shape of the triangle, from x to x' in a time St. The statistics of the 
elements of <5A are given by (f2"2"j). and Sr.; has a Gaussian distribution, characterized by (Sri) = 
and (Sri ■ Srj) = AD^SijSt, where £>b is the diffusion coefficient associated with Brownian motion. 
The scaling of the Brownian component in (|53[) implies that the increments of A and <f> due to the 
Brownian motion are given by equations (|46[) and (|48[). alternatively by stochastic incremenets 6z, 8<f) 
which can be expressed using the drift and diffusion coefficients (|52[) . Similarly, there are independent 
contributions to Sz and 8(f) from the shearing motion. These are described by stochastic increments 
with drift and diffusion coefficients given by equations ([34]) . (j35|) and (|36p . 

The stochastic terms i5A and Srt in ([53")) both act multiplicatively, either on the variables Xi, or 
upon the radius of gyration R, which is a function of these variables. This could potentially result 
in an ambiguity in defining the model depending upon whether we use the Ito or Stratonovich 
definition of the stochastic increment. We notice in this respect the ambiguity in the evolution defined 
by the Brownian term only affects the radial variable R, and not the shape variable. As pointed out 
in Sectional the same is true for the shear term. 

The fluctuations of shape are therefore described by a Fokker-Planck equation for the joint proba- 
bility density on the Kendal sphere, P(z, <j>, t) 



dP d d 2 d 2 



- [>.(*] + —i WMP] + its \D t (t)P\ (54) 



with 



v z (z) = -8D b z-8D s z 3 
D z (z) = 4(1 - z 2 )D h + Az 2 (l - z 2 )D s 

The steady state solution only depends upon the ratio of diffusion coefficients 

- 1 ■ « 

This has a steady state solution where </> is uniform on the circle and where the marginal distribution 
for z satisfies 

2z(e + z 2 )P +^-[(e + z 2 )(l - z 2 )P] = (57) 
uz 

so that 

dP 

2zP+ (e + z 2 )— = . (58) 
oz 

This has a solution 

P(z) = - _f" , r , 1 5 (59) 



2tan" 1 (l/ v / i) 



where the multiplier normalises P(z) on [—1,1]. 

To check the prediction of (|5^]l. we have studied numerically the mixed model defined by 
One of the crucial properties of the random strain is that it conserves exactly the area of the triangle: 
det(M) = 1, where M is introduced in (|18ll9p . We stress that the constraint det(M) = 1 is satisfied 
only when the evolution equation (|53[) is treated using the Stratonovich formulation One readily 
notices that the expression for the evolution operator, induced by the strain, and correct up to order 
St 1 , is simply: 

x'i = (I + 5 A + 5A 2 /2) Xl (60) 
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-j n-J I 1 1 — I — 1 1 — I — 1 1 — I — I 

10" 3 10" 2 z 10" 1 10° 

Fig. 1 The PDF of z for values of e = D b /D s , e = 3 x 10 -4 , 10" 3 , 3 x 10" 3 , 10" 2 and 3 x 10" 2 , shown as a 
full line. The fit by equation Q59p is shown as a dashed line. 



which is the expansion of the exponential of the traceless matrix SA up to second order. To impose 
strictly the constraint of area conservation, we replaced the operator in (|60[) by exp((5A), whose deter- 
minant is exactly equal to 1. and has the correct expansion, up to terms of order St 3 / 2 . Using the fact 
that: SA 2 = {8a 2 + 8b 2 )I, one readily finds that 



exp(^A) = cosh(V(<5« 2 +^ 2 ))I + S[lM j {6 f SA ■ ( 61 ) 

V(<)a 2 + 5b 2 ) 

Also, the Brownian motion term acts in ([1]) multiplicatively (because of R). As we are interested here 
only in the shape evolution, and not on the evolution of the size, R, which requires some care, we have 
projected the term RSri term in the direction transverse to the radial direction, which does not affect 
the shape dynamics, but simplifies the numerical integration. 

The data shown in Figure 1 was obtained by iterating the mapping given by (|53[) for 3 x 10 time 
steps, over 4 x 10 5 different realisations. Figure 1 shows our numerical results, which agree very well 
with the theoretical prediction. 



6 Alternative formulation using Casimir operators 

The derivation of the stationary shape distribution in sections |4][5] proceeded by setting up Langevin 
equations for the evolution of the variables z and cf> (which label positions on the Kendall sphere of tri- 
angle shapes) . This Langevin process was used to obtain a one-dimensional Fokker-Planck equation for 
z, after first concluding that the probability density is uniform in <j). An alternative approach is to derive 
a Fokker-Planck equation for the distribution of all of the variables z, </>, 6 and R which parametrise 
the relative disposition of three points in the plane, and to solve the resulting four-dimensional Fokker- 
Planck equation by separation of variables. This alternative method makes plain the invariancc proper- 
ties under the group of rotation, SO(2) and the group of dilation, and in the case of the random strain 
model, under the group of unitary matrices, SL(2) [7], at the cost of complicated algebraic manipula- 
tions, which we do not present in full detail here. These invariances are expressed via the introduction 
of the Casimir operators, L 2 , G 2 and A corresponding respectively to SO(2), SL(2) and the group of 
dilation (7j. Consistent with the convention introduced in [7], we use for the purposes of this section 
the coordinate £ = dct(w) to characterise the size of the triangle, rather than the radius of gyration 
R. In fact, triangles can be parametrized either by the set of coordinates (z, </>, x, C) coordinates, as in 
0[T], or by the set of coordinates (z, 0, x, R) used so far. Note that our notation differ from those of [7] 
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and pQ, where the coordinates (z,<f), x, C) were denoted (w, X-,9-, () (we changed the notation because 
it is conventional to use (z,4>) as the polar projection and azimuthal angle on a sphere). 



6.1 Fokker- Planck equation for random strain model 

Here we consider how to describe the random strain process using a Fokkcr-Planck equation. We 
parametrize a triangle by the two vectors, u a , a = 1, 2, introduced in equation (|3|). These vectors are 
specified by four coordinates it a ,i, and we start by describing the random strain process by a Fokker- 
Planck equation expressed in terms of these four variables. Under the action of the random strain 
process, these vectors evolve according to 

u a = a(t)u a (62) 

where a(t) is a random 2x2 matrix. The antisymmetric part of a(i) corresponds to an overall rotation. 
Since we are interested only in the shape, we ignore the overall rotation and dilation, and consider 
only traceless symmetric matrices. 

By integrating equation (|62j) with respect to time and iterating the resulting equation, in the limit 
of small values of the time St one obtains: 

««,*(«)= ««,i(0) + / dt' Oij(t'K,j(0)+ / dt' dfa^t^ak^Ua^ + OiSt^ 2 ) (63) 
Jo Jo Jo 

(here and in the remainder of this section summation is implied when indices are repeated). In equation 
(1631) . the term written as a single integral has zero mean, and a variance growing like St. It is therefore 
a diffusion term. The term written with a double integral has a mean of order St; it corresponds to the 
drift term. Higher order terms are formally smaller order terms in St. To make the connection with 
the formulation used in section [3J the stochastic increment in equations (|20j) and (|21[) may be written 



5A 



j-t+St 

J dt 1 a(i') . (64) 



We assume that the coefficients of the matrix a(<) are (5-correlated in time. Using the constraint 
that a is symmetric, traceless, and has rotationally invariant statistics, we find 

(a l:j {t)a k i{Q)) = (a 2 ) x S(t) x (5 ik Sji + SuS jk ~ (65) 

where (a 2 ) is the only parameter that characterizes the amplitude of the noise. This relation is equiva- 
lent to equations (f2~Tj) and (|22p . with (a 2 ) = 2D S . The use to the relation (|6"3")) for advancing the variables 
u a i amounts to using the physically correct Stratonovich procedure to integrate the stochastic model 
©• 

From equation (|63l) and (|65l) one can easily determine the Fokkcr-Planck equation for the random- 
strain process: 

d t P = -^—(v a ,iP) + ( " 2) „ d l ([S ik u a (0) ■ u b (0) + u a , k (0)u b ,i(0) - u a 40)u b . k (0)}p) 

ou a ,i 2 du a ^du b ^ \ J 

d (a 2 ) - 

(v a ,iP) + ^rLdP (66) 



du a .i ' 2 

where v a> i is the drift term velocity, obtained by taking the average of (|63p . and Ld is a diffusion 
operator containing second derivatives. For the drift term we find 



v a ,i = (a 2 )u a ,i 



(67) 
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The term representing the diffusion term in the Fokkcr-Planck operator has a more complex structure. 
We find it convenient to write Ld as the sum of three operators, 

Ld'P= n d l [ua(0)-ti»(0)P] (68) 

0U a ,iOU b . 3 

tl-P=o^KM- Ub Ao)P}- (70) 



6.2 Representation using Casimir operators 

Now consider how the Fokker-Planck equation may be transformed from Cartesian coordinates 
to the more convenient set (z, <f>, x, C) which separate the actions of shape, rotation, and scale trans- 
formations. As a consequence, we shall see how the operator may also be decomposed into Casimir 
operators for these fundamental group operations. 

The that drift term in the Fokker-Planck equation, dl(v l a P), makes a contribution which takes a 
very simple form, involving only the coordinate £: 

f) / BP \ 

Each of the three operators in equation (jSB"]) can be written in terms of the Casimir operators of 50(2) 
(L 2 ), of SL(2) (G 2 ) and of the dilation operator, A. 

The angular momentum operators are, up to a factor ±i, 

Lij - p =^d^- ua ' j d!- )p (72) 

and in this two-dimensional problem the corresponding Casimir operator is 

L 2 = -L^L ji . (73) 

A straightforward calculation leads to: 



where d = 2 is the space dimension, and where 



L " P = 73 5 (u ad u b>i P) - («„ ■ u b P) - (d- l)(d 2 + A)P (74) 



d 

A = u a>i - . (75) 



The operators that generate the representations of SL(2) are: 



& b = u a , i J--*£A. (76) 
ouhi 2 



The operator G ab and A commute with each other. 

The operator G 2 = G ab G ba can be easily expressed in the form: 



62 = — K<«w) - -M + d2 )( A + 2d2 ) ■ ( 77 ) 

oub,idu a j d 



Last, we express A 2 : 

d 2 



A 2 = j— ^ (u a ,iU b , r ) - (2d 2 + l)(d 2 + A)+d 4 . (78) 

du„idu b j 
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Thus, the 3 operators L d are: 



L\ = -L 2 + G 2 + -{A + d 2 ){A + 2d 2 ) - (d- l)(i + d 2 ) 

L 2 = G 2 + ^(A + d 2 )(A + 2d 2 ) 
L\ = A 2 + (2d 2 + l){A + d 2 ) -d 4 . 

With these expressions, one can simply obtain the expression for Ld'. 

L d -U = (-L 2 + 2G 2 + 2A + 8)U . 



(79) 

(80) 
(81) 

(82) 



The terms L 2 , G 2 and A 2 all have a simple expression in terms of the variables (z, cf>, xX) 12 • Specifi- 
cally, 



dx 2 



2z 2 



A 



(1-z 2 )- 
d 2 d 



z 2 ( d 2 

■I 1- 2z 



d 2 



2(1 -z 2 ) \d X 2 d(j) 2 



d 2 \ 
dxd<t>' 



(83) 
(84) 
(85) 



Thus, the operator Ld reads: 



L d -U = 4z 2 — 



dz 
1 d 2 U 

1 — z 2 d(j> 2 



dz 



d 2 2z d 2 

(l — z 2 )\d4> 2 dxdcj); 



«f + «7 



(86) 



6.3 Relation to formulation in section [3] 

Equation (j86|) enables us to express the Fokker-Planck operator in terms of the coordinates (z, 0, x, C)- 
Let us consider how to obtain the results developed in section [3] using this representation of the Fokker- 
Planck operator. Because we intend to look at solutions which are independent of the rotation angle 
X, the diffusion part of the operator is: 



n d 

L d -U = 4z 2 — 

az 



oz 



d 2 U 



(1 - z 2 ) dtp 2 



4cf + 8U. 



(87) 



Thus, the full Fokker-Planck equation, expressed as in equation (|66|) . is simply 



d t P = 



(« 2 >( 



2z< 



dz 



(l-z 2 ) 



dP 

dz 



d 2 P 



2(1 - z 2 ) d(j> 2 



As expected, this is independent of £. However, its structure appears hard to reconcile with the results 
for the random strain model in section [3] In order to make the connections clear, we first remark that 
this is an equation for a probability density P expressed as a function of the four original variables, 
u a .i, a = 1,2, i = 1,2. The transformation to the new variables C, Xi z i involves the Jacobian 



D(z,<f> >X ,0 ~ z 2 



(89) 
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The total probability is obtained by integrating P over all the u a ^. After transforming to the (z, </>, \, C) 
variables it is more convenient to deal with a probability density P'(z, (/>, x, C) which is a conserved 
quantity in these variables. In the (z, </>, x, Q variables, the normalization condition reads: 



J dtti,i diti j2 du 2 ,i du 2 ,2 P(ui,i, ui,2, "2,1, "2,2) 
= y dCdxdzd0 P(C,X,2,0)^ = J d(dxdzd(/) P'(C,x,z,t 
If we rewrite equation (|88f in terms of P' we obtain 



(90) 



az 



(I-. 2 ) 



<9z 2 P' 
dz 



d 2 P 



2(1 -Z 2 ) 



Using some elementary manipulation of the part in z of the evolution operator, one can rewrite: 

dz 2 P'~ 



d_ 

dz 



(i-z^y- 



dz 



^[, 2 d-^]-|(-2,3p) 



(91) 



(92) 



which agrees with the conclusion that the change in shape can be described by a diffusion z 2 (l — z 2 ) 
and a drift — 2z 3 . 



6.4 Fokker-Planck equation for Brownian process 



Next consider the Fokker-Planck equation for the relative Brownian diffusion of three points, and its 
relation to the derivation of Kendall's result on diffusion which as discussed in section^ After a lengthy 
calculation, this is found to take the following form 



d t P 



Db 



2C 



,d 2 P 



dC, 2 
.2 



2z 2 ^[(l-z 2 )d z P] 



Q2p Q2p 



2(1 -z 2 ) \ d(t) 2 d X 2 



2z 



d 2 p \ 

dxd<t>) 



+ 4(l-z 2 )< 



d 2 p 

dzd( 



where is the amplitude of the additive noise term: 

dtu a ,i=Va,i with (v a ,i(t)v bi j(0)} = 2D h S(t)6 ab 6ij 



(93) 



(94) 



Because we chose to consider a process that does not change the radius of the solution. This suggests 
that the evolution will be simple in the variables (z,<p,x, R)- Indeed, the evolution operator in these 
variables is simply: 



d t P = D h x 



0_ 

dz 



(l-z 2 )^-P 
az 



1 



1-z- 



■(#02 + d x 2 - 2zd<t,d x )P 



(95) 



The Jacobian of the transformation to the variables (z,(j>,x,R) is proportional to R 3 , i.e., it does not 
depend on the angular variables. The simple identity: 



d_ 

dz 



dz 



d 2 



)P] - £-{-teP) 
dz 



(96) 



immediately leads to the conclusion that the evolution of z results from a diffusion (coefficient oc 
(1 — z 2 )) and a drift (coefficient oc — 2z), as obtained in section^} 
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7 Summary 



We have analysed the simplest model which could describe the shape statistics for triplets of points 
advccted in a two-dimensional random flow. The model combines Brownian processes which describe 
a random strain with a Brownian motion of the triangle vertices. We characterised the shape using 
Kcndalll's spherical representation of the manifold of triangle shapes. The steady-state distribution is 
uniform in the azimuthal coordinate <f> but has a Lorentzian form in the polar altitude z, depending 
upon a single dimensionless parameter. 

The model considered in this work can be viewed as a simplification of the models proposed in [TJ[5] , 
which are taking into account the hierarchy of time scales in a turbulent flow. Understanding the shape 
distributions predicted by these models will require some generalization of the ideas developed in this 
work. Still, we expect that the solution discussed here provides a starting point for future analysis of 
the shape distribution of triangles, or more complicated sets of points transported by turbulent flows. 
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